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We show the equivalence of the functions G p (t) and \^{d, t)\ 2 for the "passage time" in tunneling. 
The former, obtained within the framework of the real time Feynman histories approach to the 
tunneling time problem, using the Gell-Mann and Hartle's decoherence functional, and the latter 
involving an exact analytical solution to the time-dependent Schrodinger equation for cutoff initial 
waves. 
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I. INTRODUCTION 

The tunneling time problem has remained a controver- 
sial issue after the question of how long it takes a particle 
to traverse a classically forbidden region, was raised 70 
years ago Oj. There are a number of approaches to this 
problem |2|. In this paper, an unexpected close relation- 
ship is found between a real time Feynman path integral 
approach [1 [J and the quantum shutter approach 
0)13 > which are, at first sight, unlikely to be related. 

If we use the real time Feynman path integrals 8j, 
we can define the "amplitude distribution" of tunneling 
time as the sum of e lS l h {S being the action) over the 
paths that take a specified amount of time to traverse 
the barrier region. With the amplitude distribution, we 
can deal with the interesting question whether or not a 
probability distribution is definable for tunneling time 
yj. The definability of the probability distribution de- 
pends on whether or not the amplitude distribution has 
the property of orthogonality, i.e., whether or not the 
classes of Feynman paths taking time t\ and t-i [j\ =/= T2) 
to traverse the region interfere. For rectangular barriers, 
one of the authors studied the interference quantitatively 
to conclude that (i) a probability distribution is not de- 
finable 0, but (u) the range of the values of tunneling 
time is definable In Ref. |4|, a function G(t) is intro- 
duced to analyze how different classes of Feynman paths 
(each class being characterized by the value of r) con- 
tribute to the tunneling process. The function G(t) was 
used to prove the undefinability of the probability distri- 
bution and also to estimate the range of the tunneling 
times. For typical opaque barriers, the graph of G(r) 
showed a peaked structure near the Biittiker-Landauer 
time 01 • ^ is thus clear that G(r) is an important 
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quantity for the study of the tunneling time problem. It 
is however important to understand how G(t) is related 
to the dynamics of tunneling, which is not evident at all 
from the Feynman paths construction of G(r). In the 
present paper, we will relate G(t) (to be precise, G p (r) 
as discussed below) to a time-dependent wave function. 
Now, we have to quickly add the following: In general, 
a Feynman path crosses the barrier region many times, 
so that we can define "the amount of time taken by a 
Feynman path to traverse the barrier region" in several 
ways. We can define it as the sum of the times dur- 
ing which the Feynman path is within the barrier region 
, which may be called the resident time of the Feyn- 
man path. Or, we can define it as the last time the path 
leaves the barrier region minus the first time it enters 
the region 0, which may be called the passage time of 
the Feynman path. These two different definitions at the 
level of Feynman paths would lead to physically different 
tunneling times, which we shall call the tunneling time 
of resident time type {resident time for short) and the 
tunneling time of passage time type {passage time for 
short). Reference |J concerns the resident time, while 
Refs. 0,11 and this paper concern the passage time. We 
shall attach, if necessary, subscript r to the quantities for 
the resident time (e.g., G r (r)) and subscript p for the 
passage time (e.g., G p (r)). 

Another novel approach, relevant to the tunneling time 
problem 6| , is to consider an analytic time-dependent so- 
lution to Schrodingcr's equation with the initial condition 
at t = of an incident cutoff wave, to investigate the time 
evolution of the probability density through an arbitrary 
potential barrier. This problem may be visualized as a 
gedanken experiment consisting of a shutter, situated at 
x = 0, that separates a beam of particles from a poten- 
tial barrier located in the region < x < d. At i = 
the shutter is opened, and the probability density rises 
initially from a vanishing value and evolves with time 
through x > 0. At the barrier edge x — d, the proba- 
bility density at time t, \^{x, t)\ 2 , yields the probability 
of finding the particle after a time t has elapsed. Since 
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initially there is no particle along the tunneling region, 
detecting the particle at the barrier edge at time t should 
provide a relevant time scale of the tunneling process. 
In recent work, two of the authors Q, El analyzed the 
time evolution of the probability density \^(d, t)\ 2 for a 
rectangular potential barrier using the above formalism. 
There, it was found that the probability density at the 
right barrier edge x = d, exhibits at short times a tran- 
sient structure that they named time-domain resonance. 
The maximum of the time-domain resonance, occurring 
at a time t = t p , represents the largest probability of 
finding the particle at x = d- In Ref. the above au- 
thors called the attention of the readers to the fact that 
the shape of the graph of \^(d, t)\ 2 depicted in Fig. 1 of 
that paper resembles the average shape of the graph in 
Fig. 2 of Ref. [4(, which is the graph of G r (r). Then 
they guessed that \^(d, t)\ 2 would be more related to the 
passage time rather than the resident time. In fact, in 
Ref. j5j, Yamada has studied G p (t) to find that, if r is 
simply replaced by t, the graph of G p (t) for a monochro- 
matic case is actually indistinguishable from the graph 
of \ty(d, t)/T\ 2 , where T is the transmission amplitude. 
However, there has been no explicit proof that these two 
functions are really equivalent. 

The aim of this paper is to prove that the function 
G p (r) and the probability density \^(d, t)\ 2 under the 
initial condition stated above are actually related by 
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thereby establishing a surprising relationship between the 
two approaches. As a by-product of our proof to Eq. @J, 
we present an alternative derivation, along the transmit- 
ted region, of the expression for ^/(x, t) without using the 
Laplace transform method. This derivation is the second 
purpose of the present paper. 

Section II presents a brief account of the main features 
of both approaches. Section III deals with the proof to 
Eq. (JTJ and also with a new derivation of ^(x,t). In 
Sec. IV, a numerical example is presented for a rect- 
angular potential barrier in order to exhibit the equiva- 
lence of both approaches. Concluding remarks are given 
in Sec. V. 



P is the tunneling probability defined by 



P = lim / dx^ix,^ 



(3) 



where d is the position of the right edge of the barrier. 
The decoherence functionals were formulated in general 
terms by Gell-Mann and Hartle 01 m their version of 
the consistent history approach to quantum mechanics 
[Til ITR Il6[ . The real part of D[ti;t 2 ] is a measure of 
the interference between the classes of Feynman paths 
that take different amounts of time (n and ra) to tra- 
verse the barrier region. Roughly speaking, G(r) is the 
square modulus of the sum of e lS ^ n over those paths that 
take less than time r to traverse the barrier region (to 
be precise, the result of the sum over paths is multiplied 
by the initial wave function, followed by the integrations 
over the initial and the final positions before and after 
taking the square, respectively). It is easier to deal with 
G(t) than D\ti : T2\ since G(r) is a real function of one 
variable, while D\t\; ra] is a complex function of two vari- 
ables. G(t) has the following properties: (i) G(0) = 
and (ii) G(oo) = 1. Yamada 4 : ] claimed that (a) if G(r) 
is not an increasing function of r, a probability distribu- 
tion of tunneling time is not definable, and (b) the range 
(r< , t> ) of times is an estimation of the range of values of 
tunneling time, where t< and t> are such that G(r) < e 
for v r < t< and |1 - G(r)| < e for v r > r>, where 
< e -c 1. The first claim (a) is based on the weak 
decoherence condition [TH Il6| in the consistent history 
approach. 

For a particle with wave number ko (> 0) impinging on 
the square barrier of height Vnthat extends from x = 
to x — d, G p was found to be [5j 
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ikd 
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k 2 - k 2 



. (4) 

where T(fc) = T(fc,Vo,d) is the transmission amplitude 
for the square barrier when the wave number is k, and 
T = T(k ). 



B. Quantum shutter approach 



II. THE FORMALISMS 
A. Real Time Feynman path integral approach 
In Ref. |4j, Yamada introduced G(t) by 

G(r) = i £ dn dr 2 D[ Tl -T 2 \. (2) 

In the above expression, D[t\ \ t 2 \ is the decoherence junc- 
tional for the case of tunneling time for transmission and 



A direct access to tunneling phenomena in time domain 
is to follow the time evolution of the Schrodinger's wave 
function. In Refs. 0, El , two of the authors studied the 
time-dependence of the probability density by using an 
explicit solution [6] to the time-dependent Schrodinger 
equation, with a cutoff plane wave initial condition, 



#(x,0) 



ikox 



-ik() x 



for x < 
for x > 0, 



(•5) 



impinging on a shutter placed at x = 0, just at the 
left edge of the structure that extends over the inter- 
val < x < d. The tunneling process begins with the 
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instantaneous opening of the shutter at t = 0, enabling 
the incoming wave to interact with the potential at t > 0. 
The exact solution along the transmitted region (x > d) 
reads 0, 

*(x,t) =T(k )M{x,k ;t) -T(-k )M(x,-k ;t) 

OO 

- ]T T n M(x,k n ;t). (6) 

n— — oo 

In the above expression, the quantities T(±fco) refer to 
the transmission amplitudes, the index n runs over the 
complex poles k n of T(k), which are distributed in the 
third and fourth quadrants in the complex k plane, and 
the factor T„ is defined as 



T 
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where {u„(x)} are the resonant eigenfunctions 6], which 
are the solutions to 

d 2 u n (x) \ 2 2to " 



dx 2 



u n (x) = 



(8) 



with outgoing boundary conditions, 

= -ik n u n (0), 



-r-u n (x) 
dx 



x=0 



and 



dx 



= ik n u n (d). 



(9) 



(10) 



Both the complex poles {k n } and the corresponding 
resonant eigenfunctions {u n (x)} can be calculated us- 
inga well established method, as discussed elsewhere 
013] ■ Note that from time-reversal considerations [l7j . 
the poles fc_„, seated on the third quadrant of the com- 
plex fc-plane, satisfy k— n — ~/c* and correspondingly 
U-nix) = Un(x). In Eq. jnj), the M functions are de- 
fined by 
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where q = k n , ±fc , and w(iy q ) is the complex error func- 
tion [13 with the argument y q given by 
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III. EQUIVALENCE OF BOTH APPROACHES 

A. Proof of Eq. © 

We will start from the general relationship between an 
initial wave function and the time evolved wave functions: 



dyK(x,t;v,0)9fa0), (14) 



where K (x, t] y, 0) is the propagator from (y, 0) to (x, t). 
Since our initial wave function is vanishing for x > and 
since we are interested only in the transmitted region, we 
need to know K (x, i; y, 0) only for y < and x > d, for 
which it is well-known that 



K{x,t;y,0) 



dk 
2^ 



T(k)e 
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(15) 



which follows from the cigcnfunction expansion of the 
propagator. The initial wave function can be expanded 

as 



¥fo,0) = 



dk 



(16) 



where <p(k) is the /c-space wave function. Substituting 
Eqs. (|15|) and l|l()l) into Eq. (|14fl . we can carry out the 
integration over y to have 
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For our initial wave function [Eq. (JHJ] 
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,c 1-k \k — ko + ie k + ko + ie / 
where e is an infinitesimal positive number. Thus, 

1 



(18) 



+ ie k + ko + ie / 
x T (k)e lkx - inkH / 2m \ (19) 



for x > d. 

Let us note that, since ^(x, 0) = for x > 0, 



dk 



1 



1 



k — kg + ie k + ko + ie 



T(k) e lkx = 



for x > 0, which is also apparent from the fact that the 
transmission amplitude on the complex fc-plane has sim- 
ple poles only in the lower half-plane. Owing to Eq. 12011 . 
we can rewrite Eq. I|19(l as 



1 



T{k)e ikx ( 



— kn + ie k + kn + ie 



p ihk 2 t/2m _ e ~ihko t/2m 



(21) 



Apply the following equation in Eq. 1211) 
1 _ 1 



k ± fco + ie k ± ko 



TTi6{k±ko). (22) 
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We then notice that (i) since e -^k 2 t/2 m _ e -thk 2 t/2 m = Q 
at k = kn, the contributions from the delta functions 
vanish and (ii) since (e-'" 2 '/ 2m - e -^k 2 t/2 m y( k ± k ^ 
is regular in the limit k — > Tfcoj the Cauchy principal 
value integrals can be replaced by the ordinary integrals 
(i.e., the symbol V can be removed). Consequently, we 
have for x > d 



i 

2^ 



elk 



1 



k — fcp k + ko 



T(k)e 



ikx 



^ f ^ — ihk 2 t/2m ^ — ihko 2 t/2m 



— ^— c -ihk 2 t/2m 
7T 

/oo ih(k 2 —k 2 )t/2m _ i 
dkT(k)e lkx — — g . (23) 
-oo rCo 

With this expression for ^, it is easy to see that 
|*(d,t)/T| 2 agrees with the right-hand side of Eq. (@} 
if t is replaced by t. This completes the proof of Eq. QJ. 



B. New derivation of the quantum shutter solution 

As mentioned earlier, the transmission amplitude has 
in general an infinite number of simple poles distributed 
on the lower-half of the complex fc-plane. The transmis- 
sion amplitude may be expanded in terms of its complex 
poles and corresponding residues by using a special form 
of the Mittag-LefHer theorem due to Cauchy It 



be written as 



nk)= £ 
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k k n h n J 



may 



(24) 



where r n is the residue of T(fc) at k = k n . Using Eq. 124|) 
in Eq. I|19|) . we have, for x > d, 

Z7T 



E 



dk 



o + «e 



k + kn + ie 



If we expand 
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k k n k n 
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' n \ ikx-ihk 2 t/2m 



(25) 
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1 

ko + ie k + ko + ie J \k — k n ' k n 

and use the partial fraction decompositions, we see that 
the right-hand side of Eq. I|25|) can be expressed as a sum 
of the integrals of the form of Eq. Ijllfl . The expansion 
gives four terms, which are k~^-{k ± fco + ie) -1 and 

11 1 



k ± fco + ie k — k n ±fco + k n + ie 
1 1 



k — k n k ± ko + ie J ' 



(26) 



so that Eq. (|25|l becomes, after some algebra, 

-e' 

-E 



-fco — fc'n 



_ + _^) M (a;,-fco;0 



^ fco + k n + it 
x M(x,k n ;t). 



kn it 



(27) 



In the limit e — > 0, the sums over n in the first and the 
second lines on the right-hand side of Eq. I|27(l give T{kn) 
and T(— fco), respectively [see Eq. (1241 ] . We thus obtain 

tf (a;, t) = T(k )M(x, k ; t) - T(-k Q )M (x, -k ;t) 



-2fc J2 



n— — oo 



" 2 M(x,k n ;t). (28) 

fco fcn 



Our goal here is to derive Eq. 10. In fact, Eqs. I|28|) 
and J^J are the same because of the relationship 



iu n (0)u n (d)e %k 



//,-.., (i 



(29) 



We shall prove Eq. l|29[l to conclude this section. It 
is helpful to consider the outgoing Green's function 
G + (x,x'; fc), which is the solution to 



d 2 G+(x,x';k) 
dx 2 



k,l 



2m , , 
TTV(x) 



G + (x, x'; fc) = S(x— x') 
(30) 

with the outgoing boundary conditions. We first use the 
fact that G + {x,x[\k) can be written in terms of the res- 
onant states as 20], 

n+l ' U\ u n{x)u n {x') , . 

G^{x,x;k)= 2^ 7^771 TT (0<x,x<d). 

The above expansion holds provided that the resonant 
eigenfunctions u n (x) are normalized according to the 
condition Q 



u 2 l (x)dx + i 



ul(0) + ul(d) 
2,k n 



(32) 



Next we use the fact that the transmission amplitude and 
G+(0, d; fc) are related by [2l| 



T(k) = 2ikG + (0,d;k)e 



'!■<! 



(33) 



From Eq. (|31|) . we have 

lim (fc - fc„)G+(0, d; fc) = u n (0)u n (d)/2k n , (34) 

k — >k n 

while from Eq. I|33(l together with Eq. I|24|l . we have 

lim (fc - fc„)G+(0, d; fc) = r n e lk " d /2ik n . (35) 

Equating the two results, we obtain Eq. (|29f) . 
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4.0x10 s 




The maximum of this transient structure represents the 
largest probability to find the tunneling particle at the 
barrier edge x = d. In our example, as shown in Fig. 
^ the maximum of the time-domain resonance occurs at 
t p = 5.347 fs, faster than the free passage time across the 
same distance of 10.083 nm, that is, t p /tf — 0.455. From 
t/tf = 2.0 onward the probability density approaches es- 
sentially to its asymptotic value. We have also included 
in Fig. the plot of G p (t) (dotted line) calculated from 
Eq. 0} for the same set of parameters; it is indistin- 
guishable from the previous calculation, i.e., both curves 
coincide exactly. 



V. CONCLUDING REMARKS 



FIG. 1: Plot of |*(d,i)| 2 (solid line) and G p (t) (dotted line) 
at the barrier edge x — d — 10.083 nm as a function of time in 
units of the free passage time tf. The time-dependent solution 
is normalized to the transmission coefficient \T\ 2 . 



IV. EXAMPLE 



We have found a surprising relationship between the 
real time Feynman histories approach and an analytical 
expression for the probability density for cutoff initial 
waves involving the quantum shutter setup for the "pas- 
sage time" in tunneling. This may prove to be of interest 
in the pursue of elucidating the notion of tunneling time 
through a classically forbidden region. 



To exemplify the time evolution of the probability den- 
sity we consider the set of parameters: V = 0.70 eV, 
d = 10.083 nm, E = 0.140 eV, m = 0.067 m e (m e be- 
ing the bare electron mass), inspired in semiconductor 
quantum structures |22j . In this particular example, the 
potential barrier parameters are chosen in such a way 
that k d = V /E = 5, where k = [2-mE] 1 / 2 /h. The 
opacity a of the barrier is defined as a = k'd, where 
k' = pmVo] 1 / 2 /h. In our case a — 11.18, correspond- 
ing to an opaque barrier (a 3> 1). The solid line in 
Fig. [U shows \^(d,t)\ 2 calculated with Eq. © at the 
barrier edge x = d as a function of time in units of the 
free passage time tj = md/Tik = 11.753 fs. At early 
times one sees a time-domain resonance structure Ml. 



Acknowledgments 

G. G-C. and J.V. acknowledge partial financial support 
of DGAPA-UNAM under grant No. IN101301. N. Y. is 
grateful to the Instituto de Ffsica, UNAM for their hos- 
pitality and for financial support from the Tomas Brody 
Spitz Chair. Also, N. Y. acknowledges partial finan- 
cial support of Grant-in- Aids for Scientific Research from 
the Ministry of Education, Culture, Sports, Science and 
Technology, Japan and thanks Professor H. Yamamoto 
and Professor S. Takagi for valuable discussions and to 
the Information Synergy Center at Tohoku University for 
CPU time. 



[1] L. A. MacColl, Phys. Rev. 40, 621 (1932). 

[2] See for example: E. H. Hauge and J. A. St0vneng, Rev. 
Mod. Phys. 61, 917 (1989); R. Landauer and Th. Martin, 
Rev. Mod. Phys. 66, 217 (1994); Time in Quantum Me- 
chanics, edited by J. G. Muga, R. Sala, I. L. Egusquiza 
(Springer- Verlag,' Berlin, 2002). 

[3] N. Yamada, Phys. Rev. A 54, 182 (1996). 

[4] N. Yamada, Phys. Rev. Lett. 83, 3350 (1999). 

[5] N. Yamada, (unpublished) (2002); in Meeting Abstracts 
of the Physical Society of Japan, 55, Issue 2, Part2, 
23pWD-3, 227 (2000) (in Japanese). 

[6] G. Garcfa-Calderon and A. Rubio, Phys. Rev. A 55, 3361 
(1997). 

[7] G. Garcfa-Calderon and J. Villavicencio, Phys. Rev. A 
64, 012107 (2001). 



[8] R. P. Feynman and A. R. Hibbs, Quantum Mechanics 
and Path Integrals (McGraw-Hill, New York, 1965). 

[9] Here we consider such a probability distribution of tun- 
neling times that is intrinsic to the scattering event. 
[10] M. Biittiker and R. Landauer, Phys. Rev. Lett 49, 1739 
(1982). 

[11] D. G. Sokolovskii and L. M. Baskin, Sov. Phys. Tech. 
Phys. 30, 1076 (1985). 

[12] L. S. Schulman and R. W. Ziolkowski, in Proceedings of 
Third International Conference on Path Integrals from 
me V to Me V, edited by V. Sa-yakanit et al. (World Sci- 
entific, Singapore, 1989), p. 253. 

[13] G. Garcfa-Calderon and J. Villavicencio, Preprint arXiv: 
quant-ph/0008014 (2000). 

[14] M. Gell-Mann and J. B. Hartle, in Proceedings of the 3rd 



6 



International Symposium on the Foundations of Quan- 
tum Mechanics in the Light of New Technology, edited 
by S. Kobayashi et al. (Physical Society of Japan, Tokyo, 
1990), p. 321. 

[15] M. Gell-Mann and J. B. Hartle, in Proceedings of the 25th 
International Conference on High Energy Physics, edited 
by K. K. Phua and Y. Yamaguchi (World Scientific, 
Singapore, 1991), p. 1303. 

[16] J. B. Hartle, Phys. Rev. D 44, 3173 (1991); M. Gell-Mann 
and J. B. Hartle, Phys. Rev. D 47, 3345 (1993). 

[17] J. Humblet and L. Rosenfeld, Nucl. Phys. 26, 529 (1961). 

[18] V. N. Faddeyeva and N. M. Terent'ev, Tables of values 
of the function w(z) = e~ z (1 + J Q Z e* dt) for com- 
plex argument, translated from the Russian by D. G. Fry 
and B. A. Hons (Pergamon, London, 1961); Handbook 



of Mathematical Functions, edited by M. Abramowitz 
and I. A. Stegun (Dover, New York, 1965), p. 297; G. 
P. M. Poppe and C. M. J. Wijers, ACM Transactions 
on Mathematical Software, 16, 38 (1990); 16, 47 (1990). 
The above references discuss also numerical methods to 
evaluate the complex error function w. 

[19] E. C. Tichmarsh, The Theory of Functions (Oxford Uni- 
versity Press, London, 1939), 2nd ed., p. 110. 

[20] G. Garcfa-Calderon (unpublished). A similar expansion 
in 3 dimensions is presented in the paper, G. Garcia 
Calderon and A. Rubio, Nucl. Phys. A 458, 560 (1986). 

[21] G. Garcfa-Calderon, Solid State Commun. 62, 441(1987). 

[22] D. K. Ferry and S. M. Goodnick, Transport in Nanos- 
tructures (Cambridge University Press, United Kingdom, 
1997), pp. 91-201. 



